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Abstract 

We derive the exact formula for thermal-equilibrium spacing distribution of one-dimensional particle gas with re- 
pulsive potential V(r) = r~ a (a > 0) depending on the distance r between the neighboring particles. The calculated 
distribution (for a = 1) is successfully compared with the highway-traffic clearance distributions, which provides a de- 
tailed view of changes in microscopical structure of traffic sample depending on traffic density. In addition to that, the 
observed correspondence is a strong support of studies applying the equilibrium statistical physics to traffic modelling. 

PACS numbers: 05.70-a, 05.20-y, 45.70 Vn 
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Investigation of one-dimensional particle ensembles 
seems actually to be very useful for understanding of the 
complex system called vehicular traffic. Beside the favorite 
cellular automata, in the recent time new trend appears 
in the traffic modelling. Application of the equilibrium 
statistical-physics to traffic ensembles (queuing systems 
of spatially interacting particles) has been more times dis- 
cussed (see for example Ref. |2j and other references in) 
and successfully demonstrated in the articles (5), (9), and 
Q. It has been manifested in Ref. |9| that thermodynamics 
approach can be applied to such a many-particle driven 
system as traffic flow, based on a microscopic description, 
in analogy to equilibrium physical systems. Besides, it has 
been demonstrated in Ref. J7| that the relevant statisti- 
cal distributions obtained from the local thermodynamical 
model are in accord with those from traffic data measured 
on real freeways by the induction-loop detectors. It opens 
up an opportunity for finding the analytical form of proba- 
bility density for clear distance among the cars (i.e. clear- 
ance distribution in traffic terminology) moving in traffic 
stream. Up to the present time the question about analyti- 
cal form of clearance distribution has only been a subject 
of speculations (see the Ref. (6), 13) - never successfully 
answered. 

We aim to use one-dimensional thermodynamical parti- 
cle gas to the prediction of microscopical structure in traffic 
flows and consequently make a comparison to the relevant 
traffic data distributions. Justification for the approach de- 
scribed can be found in Ref. Q and where it is proved 
that equilibrium solution of certain family of the parti- 
cle gases (exposed to the heat bath with the temperature 
T > 0) is a good approximation for steady-state solution 
of driven many-particle systems with asymmetrical inter- 



actions. Since the vehicular traffic is a dissipative system 
of active elements (moving far from equilibrium), it is evi- 
dent (see Ref. |9| and 0) that thermodynamics-balance 
approach can be used on a mesoscopic level only (i.e. for 
small-sized samples of N vehicles), where the traffic den- 
sity fluctuates around the constant value and distances and 
velocities are mutually uncorrected. Then the stationary 
solution of relevant Fokker-Planck equation practically co- 
incides with thermal-balance probability density of certain 
statistical gas (see Ref. Q). Moreover, the possibility for 
using the thermodynamical approach is supported by the 
fact that distribution of velocities in traffic sample fully 
corresponds to the Gaussian distribution (in each traffic- 
density interval) - see Ref. 13, J4), and Q. Overall, it is 
sufficiently justified that traffic systems can be locally (i.e. 
on a mesoscopic level) described by instruments of equi- 
librium statistical physics. 

Thus, consider N identical particles (vehicles) on the 
circle of the circumference L = N. Let &j (i = 1 . . . N) 
denote the circular position of the ith particle. Put xn+i = 
x\ + 2ir, for convenience. Now we introduce the short- 
ranged potential energy 
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where Vfa) corresponds to the repulsive two-body poten- 
tial depending on the distance r.; = — Xi\^- between 
the neighboring particles only. Nearest-neighbor interac- 
tion is chosen with the respect to the realistic behavior of 
car-driver in traffic sample (see Ref. J71)- Besides, the po- 
tential V(r) has to be defined so that lim r ^o+ V( r ) = 00 
which prevents particles passing through each other. The 
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hamiltonian of the described ensemble reads as 
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with the ith particle velocity i>j and the positive constant C. 
Note that v represents the mean velocity in the ensemble. 
Then, the appropriate partition functionj 
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leads us (after 2N — 1 integrations) to the simple assertion 
that velocity v of particles is Gaussian distributed, i.e. 



P(v) 
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is the corresponding probability density. 

Of larger interest, however, is the spacing distribution 
Pp(r). To calculate the exact form of Pp{r) one can re- 
strict the partition function ([T]i by N velocity-integrations 
to the reduced form 



Z N {L) 
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where (3 = Ccr~ 2 is the inverse temperature (dimension- 
less) of the heat bath. Denoting /(r) = e~P v ( r > the previ- 
ous expression changes to 

Z N {L) = [ <J|L-5>III /( r «) dr i • ■ ■ drN - 

Applying the Laplace transformation (see the Ref. UJ for 
details) one can obtain 
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(p)= / Z N (L)e-P L dL = 
Jo 



N 



f(r)e'P r dr\ = [g(p)] 



Then the partition function (in the large N limit) can be 
computed with the help of Laplace inversion 
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g N {p) e Lp dp. 



Its value is well estimated by the approximation in the sad- 
dle point B which is determined using the equation 

g(B)d P [ ' N' 



Thus, 



Z N (L) « [g(B)r e 
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Hence the probability density for spacing n between the 
particles 1 and 2 can be then reduced to the form 

ut \ Z N -i(L-i\) 

P(ri) = ^—<T\ ^( ri )- 

Z N (L) 

Supposing N 3> 1 and using equation (O we obtain 

^) = ^y/(n)e--S 

which leads (after applying the same procedure for every 
pair of successive particles) to the distribution function for 
spacing r between arbitrary couple of neighboring particles 



Pp{r) = Ae-^^ e~ Br (r > 0). 



(3) 



Note that constant A assures the normalization 
J °° Pp(r) dr = 1. Furthermore, returning to the origi- 
nal choice L = N, the mean spacing is 



e / rPp(r) dr = 1. 
lo 



(4) 



Two above conditions can be understood as equation sys- 
tem for unknown normalization constants A, B. 

Let us to proceed to the special variants of the gas stud- 
ied. Firstly, we draw our attention to the Coulomb gas with 
the logarithmic potential 



V(r) := - ln(r) 



(r > 0). 



Such a gas (usually called Dyson's gas, for example in Ref. 
UT]) is frequently used in the many branches of physics 
(including the traffic research in Ref. [8]) and the corre- 
sponding spacing distribution reads as (see Ref. UJ) 



P R ( r ) = r /3 e -(/m)r 

m ' r(/3 + i) 



where T(£) is gamma function. Of larger physical interest, 
as demonstrated in Ref. (5) and Q, seem actually to be the 
potentials 



V a (r) 



(r > 0), 



for a > 0. The aim of the following computational proce 
dure is to normalize the distribution 
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Consider now the favorable choice a = 1, for which the 
normalization integrals are exactly expressed as 
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dr 



V ~~B 
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where JC\ is the Mac-Donald's function (modified Bessel's 
function of the second kind) of order A, having for A = 1 
and A = 2 an approximate expression 
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1 <T is the constant representing a statistical variance 
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Applying the equations (O and (jT]) to the normalization in- 
tegrals one can determine the exact values of the constants 
A and B. Both of them can be, after applying Taylor's ex- 
pansion procedure, very well estimated by the approxima- 
tions 

/?+^4 — . (8) 
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and 
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Figure 1. Normalization constant B depending on the inverse 
temperature p. 

Squares represent the exact value of B obtained from numerical 
computations. The dashed and solid curves display the large P 
approximation tlOt and full approximation (8), respectively. The 
behavior close to the origin is magnified on the inset. 

Finally, we investigate the distribution (0 for general 
a > 0. Although in this case the normalization integrals are 
not trivially solvable, the scaling <j4j leads us to the simple 
approximate formula 

a 
2 

The large f3 estimation r~ a m 1 — a + or -1 , which holds 
true for values r around mean distance r» 1, provides the 
asymptotical formula for normalization constant A : 



B^a/3 + l + - 
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Figure 2. Relative deviation in the approximate value of the 
normalization constant B as function of parameter p. 

We display the deviation i\3\ between the numerical value B ex and the 

value B est obtained from the large P approximation j 10| . The plus 
signs, squares, triangles, circles and stars correspond to the parameters 
a = 1, 2, 3, 4, 5, respectively. The tails of the curves are magnified on 
the inset. 



For practical applications it seems to be useful to detect 
the critical inverse temperature (3 cr it under which the rela- 
tive deviation between the exact (ex) and estimated (est) 
values of the constant A (or B) 



'log(A) 



\\og(A ex ) -log(A est )\ 



log(Aex) 
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are larger then the fixed acceptable deviation S. For these 
purposes we plot the functional dependence 5b = 5b{P) 
and r5i og (A) = <5io g (A)(/?) in the Fig. 2 and Fig. 3, respec- 
tively. We note that the exact values A ex , B ex were deter- 
mined with the help of the numerical computations. 




Figure 3. Relative deviation in the approximate value of the 
normalization constant A as function of parameter 8. 

Plotted is the deviation H2\ between the numerically computed value 
A ex and the estimated value Jilt . The symbols used here are consistent 
with the symbols in the Fig. 2. 

Considering now the tested choice for car-car potential 
V(r) = r^ 1 we intend to compare the equilibrium distri- 
bution 



Pp(r) = Ae ~e 



-Br 



(14) 



with the relevant distributions of single-vehicle data mea- 
sured continuously during approximately 140 days on 
the Dutch two-lane freeway A9. Macroscopic traffic den- 
sity q was calculated for samples of N = 50 subse- 
quent cars passing a detector. For the purposes described 
above we divide the region of the measured densities q £ 
[0,85 veh/km/lane] into 85 equidistant subintervals and 
separately analyze the data from each one of them. The 
sketched procedure prevents the undesired mixing of the 
states with the different inverse temperature (3, i.e. with the 
different density. Bumper-to-bumper distance r, among the 
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succeeding cars (ith and (i — l)th) is calculated (after elim- 
inating car-truck, truck-car, and truck-truck gaps) via stan- 
dard formula 

n = Vi{u -U-i), 

by means of netto time-headway ij — and velocity Vi of 
ith car (both directly measured with induction-loop detec- 
tor) supposing that velocity Vi remains constant between 
the times tj,tj_i when ith car and the previous one are 
passing a measure point. Such a condition could be ques- 
tionable, especially in the region of small densities where 
the temporal gaps are too large. However, the influence of a 
possible error is of marginal importance, as apparent from 
the fact, that distribution function plotted for small-density 
data do not show any visible deviation from Poisson behav- 
ior expected for independent events (see Fig. 4 and relation 

We note that mean distance among the cars is re-scaled 
to one in all density-regions. The thorough statistical analy- 
sis of the traffic data leads afterwards to the excellent agree- 
ment between clearance distribution computed from traffic 
data and formula (TBI for fitted value of inverse temperature 
(3 fa (see Fig. 4). We have obtained the fit parameter /3fa 
by a least-square method, i.e. minimizing the error function 
X 2 ■ The deviation \ 2 between the theoretical and empirical 
clearance distributions is plotted in Fig. 5 (low part). 
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Scaled clearance r 
Figure 4. Probability density P(r) for scaled spacing r between 
successive cars in traffic flow. 

Histograms represent the clearance distributions computed for traffic data 
from indicated density region (in vch/km/lanc). Note that mean 
distance among the cars is re-scaled to one in all density-regions. The 
curves represent the predictions of statistical model H4\ for fitted value 
of inverse temperature (3 fa- The respective values Pfu were carefully 
analyzed and consecutively visualized in the following Fig. 5 (top part). 

Dimensionless inverse temperature (3 of the traffic 
sample, representing a quantitative description of mental 
strain under which the car-drivers are in a given situa- 
tion, shows non-trivial dependence on traffic density g (as 
visible in Fig. 5 - top part). For free flow states (g < 
20 veh/km/lane) one can recognize a rise in temperature 
having a linear behavior (up to 10 veh/km/lane) and vis- 
ible plateau above. In the intermediate region (between 20 
and 50 veh/km/lane), where free traffic converts to the 
congested traffic, we detect a sharp increase in the first half. 



Such a behavior can be simply elucidated by the fact that 
the drivers, moving quite fast in a relatively dense traffic 
flow, are under a substantial psychological pressure, which 
finally results (for densities g E [35, 50] veh/km/lane) 
in the transition to the congested flows a therefore in the 
drop in inverse temperature. In synchronized traffic regime 
(g > 50 veh/km/lane) the drivers vigilance rapidly grows 
up which culminates by the traffic-jam formation. 
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Figure 5. Inverse temperature /3f it and deviation x 2 as a function of 
traffic density g. 

Squares stand for values of fit parameter (3 fit, for which the function 
U4) coincides with clearance distribution of traffic data. The curve 
represents a polynomial fit of relevant data. Bars from lower part 
correspond to the sums of squared deviations between the empirical and 
the theoretical netto distance distributions for (3 fa- We note that the 
value of normalization constant B was determined via formula {8}, 
because the values of relevant inverse temperatures lie in interval [0, 2], 
where the linear approximation )10t is less suitable (see deviations in 
Fig. 2). Second normalization constant A was calculated by means of 
equation |9}. 

For completeness we have compared real-road clear- 
ance distributions with probability density © specified for 
power-law potentials V(r) = r~ a , where a = 2,3,4,5, 
respectively. Similar analysis was already introduced in ar- 
ticle 0. As results from the careful analysis of statistical 
deviations \ 2 , relevant distributions (O - obtained by ap- 
plying a least-square method - predicate the inter-vehicle 
gaps substantially worse than those calculated for a = 1 
(follow Fig. 6). We emphasize that small deviations x 2 (in- 
dicating a good agreement) detected near the origin in Fig. 
6 are caused by the fact that for low traffic-densities the in- 
teractions among the cars are vanishing and inter-vehicle 
gaps are therefore practically independent. The relevant 
distributions in this case came close to the Poisson distri- 
bution (see upper left-hand subplot of Fig. 4). However, the 
Poisson distribution can be obtained as a limit of distribu- 
tion (01, i.e. 



lim Pp{r) = e 



(15) 



for arbitrary function V{r). For that reason it is impossible 
to detect the interaction potential in traffic sample using 
the low density data only. The comparison of various po- 
tentials V(r) = r~ a (including the logarithmical potential 
V(r) = — ln(r)) brings finally the message that interac- 
tion among the vehicles in traffic stream can be very well 
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estimated by the short-ranged two-body power-law poten- 
tial V(r) = r -1 . Predicted inter- vehicle-gap distributions 
correspond in this case to the computed probability density 
in all density intervals. 
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Figure 6. Deviation \ 2 between clearance distribution of cars 
moving in traffic stream and normalized distribution {D evaluated 

for a = 2, 3, 4, 5, respectively. 

Grey bars represent the sums of squared deviations depending on traffic 
density for a = 2, 3, 4, 5, respectively. Dark bars display the relevant 
deviations for a = 1 previously plotted in Fig. 5 (lower part). 

We append that another suitable quantity for compar- 
ison with single-vehicle data is a time-clearance distribu- 
tion as well. Furthermore, time gaps among the succeeding 
cars are directly measurable and therefore not burden with 
errors caused by computation approximations. Determina- 
tion of exact form for time-headway probability density of 
above thermodynamical traffic gas and relevant compari- 
son with highway traffic data will be incorporated in the 
continuing work. Nevertheless, several analytical forms for 
corresponding distribution have already been discussed in 
article (5| and also in book lHOl . 



To conclude we have found the analytical form of 
the thermal-equilibrium spacing distribution for one- 
dimensional traffic gas which neighboring particles are 
repulsed by the two-body potential V = r~ a , where r 
is their mutual distance. The values of two normalization 
constants were successfully estimated by the convenient 
approximations. The calculated distribution (for a = 1) 
with one free parameter (inverse temperature (3) has been 
compared to the distance clearance distribution of the free- 
way traffic samples with the excellent outcome. It was 
demonstrated that inverse temperature of traffic sample 
non-trivially depends on the traffic density. The obtained 
agreement between experimental and calculated distribu- 
tions confirms a convenience of the traffic potential used 
for description of local traffic interactions. Presented ar- 
ticle crowns the quest for the mathematical formula for 
probability density of mutual clearances among the cars in 
traffic stream and supports the possibility for applying the 
equilibrium statistical physics to the traffic modelling. 
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